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' Corrections from Breit interaction to basic properties of atomic ^^^Cs are determined in the 

I framework of third-order relativistic many-body perturbation theory. The corrections to energies, 

■ hyperfine-structure constants, off-diagonal hyperfine 65 — 7S amplitude, and electric-dipole matrix 

elements are tabulated. It is demonstrated that the Breit corrections to correlations are comparable 
i—j", to the Breit corrections at the Dirac-Hartree-Fock level. Modification of the parity-nonconserving 

(PNC) 65 — 7S amplitude due to Breit interaction is also evaluated; the resulting weak charge 
of ^''■^Cs shows no significant deviation from the prediction of the standard model of elementary 
' particles. The neutron skin correction to the PNC amplitude is also estimated to be -0.2% with 

an error bound of 30% based on the analysis of recent experiments with antiprotonic atoms. The 
present work supplements publication [A. Derevianko, Phys. Rev. Lett. 85, 1618 (2000)] with a 
discussion of the formalism and provides additional numerical results and updated discussion of 
parity violation. 
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■ I. INTRODUCTION 

! The most accurate to date measurement of parity-nonconservation (PNC) jl], || in atoms has been carried out 
by Wieman and co-workers using ^^"^Cs ^. The observed weak charge of the nucleus, Qw, is determined as a 
combination of the experimental PNC amplitude -EpNC of the 65'i/2 — 75*1/2 transition and a theoretical atomic- 
structure parameter -EpNc/Qw- Such determined Qw provides powerful constraints on possible extensions to the 
standard model (SM) of elementary particles. The achieved precision in experiments ^ is 0.35%; however the 
required atomic structure parameter has been calculated only with an accuracy of about 1% limiting the 

accuracy of determination of the weak charge. Presently it is understood |Q, ^, ^ that a detailed account of the Breit 
corrections to basic atomic properties is required to reach the next level of precision in ab initio relativistic calculation 
• of PNC amplitudes. In particular, the Breit correction to the 65*1/2 ~ 75*1/2 PNC amplitude in ^^'^Cs accountsj'or a 

■ dominant part of the deviation |^ of determined weak charge from the prediction of the standard model 
QQ The purpose of this paper is to provide a detailed discussion of the formalism employed in Ref. ||^, ^ ^ 

tabulate additional numerical results. Since the publication of Ref. |l^ several calculations of the Breit correction 
to properties of cesium atom have been carried out |[ |ll| and a comparison between different approaches is also 
' presented here. We also calculate a value for the neutron "skin" correction to the PNC amplitude based on the 
c/3 I analysis of experiments with antiprotonic atoms [|2| . 
. ^ ' The major difference between the present analysis and earlier works on the Breit interaction in multi-electron 
c/2 ' atoms |lj, |l5|, |l^, |l^, |lj, |l^ is the systematic treatment of correlation effects, i.e. contributions beyond self- 
consistent Breit-Coulomb-Hartree-Fock (BCHF) formulation ^ These correlation effects are estimated here 
in the framework of relativistic many-body perturbation theory. It is demonstrated that these additional contributions 
are comparable to the lowest-order BCHF corrections for almost all considered atomic properties. 

The paper is organized as follows. In Secti on ^ we describe the employed many-body formalism. Numerical 



lOl and to 



results are tabulated and discussed in Section HI. We consider Breit corrections to energies, hyperfine-structure 



X 

Cq ' parity non-conserving 65 — 75* amplitude is also evaluated in Section [V 



constants, off-diagonal hyperfine 65* — 75* amplitude, and electric-dipole matrix elements. The Breit correction to 



II. METHOD 

The Breit interaction [po| |2l| , is a two-particle interaction caused by an exchange of transverse photons between 
atomic electrons. Qualitatively, it describes a magnetic interaction between electrons (so-called Gaunt interaction) 
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and retardation effect. Its low-frequency form in the Coulomb gauge, employed here, is given by-"^ 

i<j -I 

where a are Dirac matrices and is a distance between electrons. The omitted frequenc y d ependence constitutes 
as little as 1-2% of the total Breit correction to the energies for the atomic ground states |17|. The present goal of 
ah initio relativistic calculations in Cs is to reach an overall accuracy of 0.1%. Since the Breit corrections to basic 
properties of Cs are below 1%, we ignore the frequency dependence in the present analysis. It is worth noting that 
a consistent inclusion of the frequency dependence in the Breit interaction would require simultaneous treatment of 
QED self-energy correction JTsf. 



A. Many-body perturbation theory and Breit interaction 

The many-body Hamiltonian of an atomic system can be generally represented as 

H = H^ + T = Y,h^(^) + \Y.t{i,j), (2) 

i ij 

where 

ha{i)^c{a,-p,)+Pc^ + Vnuc{i) (3) 

is a one-particle Dirac Hamiltonian for an electron including Coulomb interaction with the nucleus Viiuc(*) and t{i,j) 
represents two-particle interactions. To effectively minimize the perturbing two-particle interactions one introduces a 
potential U (i) and rewrites the Hamiltonian as 



i ij i J 

For atoms with one valence electron v outside a closed-shell core a many-body wavefunction {"^v) hi the independent- 
particle approximation is a Slater determinant constructed from core and valence single-particle orbitals (j^i. These 
orbitals satisfy one-particle Dirac equation 

{ho + U) = . (5) 

The potential U is usually chosen to be spherically-symmetric and label i is a list of conventional quantum numbers 
{rii, ji,li,mi} for bound states, with rii replaced by Si for continuum. With the complete set of single-particle states 
(f>i, the Hamiltonian, Eq.(^), can be recast into the second-quantized form 

i ij ijkl 

Only certain combinations of positive- and negative-energy solutions of the Dirac equation ( |5[) a re retained in rela- 
tivistic many-body Hamiltonian (no-pair approximation p3|). The reader is directed to Rcf. ]24|, p5| and references 
therein for a detailed discussion of the problem of negative-energy states. 

We follow a convention of Ref. and label core orbitals as a, 6..., excited (virtual) orbitals as m, n, and 
valence orbitals as v,w. Indexes i,j,k,l range over both core and virtual orbitals. In this notation the lowest-order 
wavefunction is |^'«)'-°'' = |Ocorc), where quasi-vacuum state |Ocorc) = (Ilaecorc '^I) |0) represents a closed-shell 
atomic core. Introducing normal form of operator products, defined with respect to |Ocore) one can rewrite a 

two-particle operator T as a sum of zero-, one-, and two-body contributions [E6[ 



Unless specified otherwise, atomic units h = \e\ = me = 1 are used throughout the paper. 
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ijkl 

with tij = '^g^{tiaja ~ Uaaj)- In this notation the Hamiltonian reads 

l:^J{U,-U,,}■.a\a,■.+ (7) 
1 Y.ijki tijki ■ alo-l^iak ■■ ■ 

Zero-body contribution to the total Hamihonian H has been discarded since it does not affect the properties of 
valence states. It is worth emphasizing that the Breit and Coulomb interactions are of course two-particle operators; 
reference to zero-, one-, and two-body parts arises due to the separation into the normal forms of operator products 
and is just a matter of convenience. 

In the case at hand, the two-particle interaction 

T = C + B 

is a sum of the instantaneous Coulomb interaction C — X]i<j T~ ^^^^ the Breit interaction _B, Eq. (^. Corresponding 
two-particle matrix elements are designated as Cijki and bijki- The Coulomb interaction dominates and we distinguish 
two possibilities in defining the effective potential U in Eq. (ph: traditional Coulomb-Hartree-Fock (CHF) potential 
U'^^^ and Breit-Coulomb-Hartree-Fock (BCHF) potential U^'^^^ , where the Breit and Coulomb interactions are 
treated on the same footing. To differentiate between the two resulting eigensystems of Eq. (|^) we will add bar to the 
quantities pertaining to the BCHF case, e.g. £i,ai,a|. 

The conventional Coulomb-Hartree-Fock (CHF) equation reads 

{ho + f/^"^) 0, = £.0. , (8) 

[/CHF being mean-field Hartree-Fock potential; this potential contains direct and exchange Coulomb interactions of 
electron i with core electrons. A set of CHF equations is solved self-consistently for core orbitals; valence wavefunctions 
and energies are determined subsequently by "freezing" the core orbitals. The Breit-Coulomb-Hartree-Fock (BCHF) 
approximation constitutes introduction of the Breit interaction simultaneously with the Coulomb interaction into the 
above CHF equation 

{ho + t/^^"^) = £.<^. . (9) 
Compared to the CHF equations, energies, wave-functions, and the Hartree-Fock potential are modified. We discuss 



a relation between CHF and BCHF methods and the associated relaxation effect in Section II B. 

To simplify the second-quantized Hamiltonian, Eq. (|^), we use the fact that matrix elements of the Hartree-Fock 
potentials are {(j)i\U'~^^^ = Cij and {(l)i\U^'~^^^\(t>j) = Cij -\-bij. In the Coulomb-Hartree-Fock case the Hamiltonian 
reduces to a sum of the conventional Coulomb Hamiltonian 

H'c ^^Ei : ajui : +i ^ Cijki : aja^aiak : (10) 

ijkl 



and the Breit correction 



SbH'c = ^fey : ajuj : +^^6ijfei : ala'^^aiUk : • (11) 

ij ijkl 

In the case of equivalent treatment of the Breit and Coulomb interactions (BCHF case) the corresponding Hamiltonian 
is less complicated 

H'cj^B = ^ : aja^ : ^(cyfez -I- bijki) ■ a|a]aiafc : , (12) 

i ijkl 

since the effective one-body Breit term in Eq. ( |ll] ) has been "transformed away" by a proper choice of one-particle 
states. 
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FIG. 1: Sample second- and third-order Brueckner-Goldstone diagrams representing many-body contributions to matrix 
elements. Diagrams (a) and (b) arise in the random-phase approximation and (c) is the Brueckner-orbital correction. 
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FIG. 2: Representative second-order contribution to an energy of valence electron. 



Of course, finding a solution of the Shrodinger equation even with the traditional many-body Coulomb Hamiltonian, 
Eq. (10), is a nontrivial problem. Many-body perturbation theory has proven to be very successful in treating 
contributions beyond the Hartree-Fock level. In particular, ab initio relativistic many-body calculations for alkali- 
metal atoms have been performed by Notre Dame and Novosibirsk (Sydney) groups. These and other calculations 
have been reviewed recently in Ref. |2^. An accurate description of the correlations (i.e. contributions beyond 
Hartree-Fock value) plays a crucial role in high-precision calculations. One of the most striking examples of the 
importance of correlations in ^'^^Cs is the magnetic-dipole hyperfine-structure (HFS) constant A of 5Dc,/2 level. Here 
the Coulomb-Hartree-Fock value, -1-7.47 MHz, has a sign opposite to that of experimental value -21.24(5) MHz from 
Ref. ||2^. The dominant correlation corrections to matrix elements arise because of core-shielding of externally 
applied fields (e.g. nuclear fields for HFS constants) and an additional attraction of a valence electron by an induced 
dipole moment of the core p9| . The former effect is described by contributions beginning at second order (random- 
phase approximation (RPA)) and the latter in third order (Brueckner corrections) of many-body perturbation theory. 
Representative many-body diagrams are shown in Fig. ^ Qualitatively, the Breit correction to a certain Coulomb 
diagram is proportional to the value of the Coulomb diagram. Therefore, in addition to lowest-order corrections we 
consider the Breit contributions to the dominant RPA and Brueckner diagrams. It will be demonstrated that these 
correlated Breit corrections in many cases are comparable to the lowest-order ones. In the BCHF basis the correlated 
Breit correction to valence energies appears in the second-order; a sample diagram is drawn in Fig. ^. The corrections 
to this class of diagrams is comparable to the modification at the Hartree-Fock level. 

In the present work we employ many-body perturbation theory (MBPT). Explicit expressions for contributions to 
matrix elements up to the third order were tabulated by Blundell et al. These authors provide formulas for 

a general perturbing potential with one- and two-body parts. In the second and third orders there are 31 distinct 
diagrams involving one-body part of the perturbation and 28 diagrams containing only two-body part. Certainly 
calculations of the Breit corrections are less complicated in the Breit-Coulomb-Hartree-Fock basis, where the one- 
body perturbation is absent. Another advantage of the BCHF basis is an automatic inclusion of important relaxation 
effect discussed below. An adequate account for the relaxation effect correction in the CHF basis would have required 
fifth-order calculations for matrix elements. 

The generalization of MBPT expressions to a simultaneous treatment of Coulomb and Breit interactions is straight- 
forward: Coulomb interaction lines are replaced by a sum of Coulomb and Breit interactions and particle (hole) lines 
by Breit-Coulomb-Hartree-Fock states (see Figs. |] and p.) Together with the corrections linear in the Breit interaction 
such approach introduces terms nonlinear in the Breit interaction. Strictly speaking these nonlinear terms have no 
meaningful theoretical basis and therefore have to be omitted. However, the Breit contribution to atomic properties 
is relatively small and the much smaller terms nonlinear in the Breit interactions can be neglected at the present level 
of accuracy. 
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FIG. 3: Comparision of Hartree-Fock and first-order relative one-body Breit corrections to (a) energies of valence states, (b) 
hyperfine-structure (HPS) constants A, and (c) electric-dipole transition amplitudes. The dotted and striped bars represent 
first-order and Hartree-Fock corrections, respectively. The relative corrections to the energies and HFS constants are defined 
with respect to experimental values, and electric-dipole amplitudes with respect to Coulomb-Hartree-Fock values. 



B. Relaxation effect 



In the CHF basis the first-order corrections to valence energies e„ and matrix elements Z^v due to the one-body 
part of the Breit interaction are given by 

Sev = byy , (13) 

OZtut, h . 

Similar one-body Breit corrections can be calculated as differences between lowest-order values found in the Breit- 
Coulomb- and Coulomb-Hartree-Fock approximations 

feHF ^ e„ _ £^ , (14) 

In Fig. ^ we present a comparison of the lowest-order one-body Breit corrections to valence energies, hyperfine- 
structure (HFS) constants A and electric-dipole transition amplitudes. The dotted and striped bars represent first- 
order, Eq. (^), and Hartree-Fock corrections, Eq. ([l^), respectively. There is a striking discrepancy between the two 
corrections for all these quantities. For example, the first-order correction to the 6s HFS constant is -0.4%, while at 
the Hartree-Fock level the Breit correction almost vanishes. 

These large discrepancies are explained by a "relaxation" effect, i.e. modification of the Hartree-Fock potential 
through adjustment of core orbitals jl^. To illustrate this effect we rewrite the BCHF equation, Eq. (^), as 

(/lo + C/°"^ + AC/) 0, = £-0, , (15) 

where the perturbing potential is AC/ = [/^chf _ [/CHF^ Further 0j — (pj + Xj, where Xj is a correction to a CFH 
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FIG. 4: Diagrammatic representation of Eq. (^. Here horizontal solid (wavy) lines represent the Coulomb (Breit) interactions. 
Double horizontal lines are expansion coefBcients ^ and little "stumps" indicate places where particle or hole lines are to be 
attached. 



wavefunction (jjj due to the Breit interaction. In the lowest order these corrections can be expressed as 



To the first order in the Breit interaction 



At/ {l)~Y.j (2) b (1, 2) 0, (2) dT2 + 



+ Ea/xi(2)c(l,2)0„(2)dr2+ (17) 
+ Ea / (2) c (1, 2) xa (2) dT2 + exchange. 

Here the first term is an explicit Breit contribution, while in the second and the third terms the Breit interaction 
enters implicitly through corrections to core orbitals. Only the first term (and its exchange form) are included in the 
first-order correction, Eq. (p^. As demonstrated by Lindroth et al. for Breit corrections to the energy levels of 
Hg, the residual "relaxation" terms are large and substantially modify the first-order corrections. Similar observation 
has been made by Johnson et al. in calculations of Breit corrections to energies of sodium-like ions. Independent 
to the present analysis ( partially published in Ref. Q) the relaxation effect in Cs has been recently discussed by 



Kozlov et al. |32 



At this point it is clear that the inclusion of Breit interaction in the Coulomb-Hartree-Fock equations (i) greatly 
simplifies many-body perturbation expansions, and (ii) automatically accounts for the significant relaxation effects. 
In other words, compared to the traditional Coulomb-Hartree-Fock formulation, the transformation to the Breit-CHF 
basis sums many-body diagrams involving the effective one-body Breit interaction to all orders of perturbation theory. 



C. Construction of Breit-Coulomb Hartree-Fock basis 



Several methods can be devised for constructing the Breit-CHF single-particle basis. For example, one can determine 
the Breit corrections to wavefunctions by substituting Eq. ([l7| ) into Eq. (|l6|). It is convenient to express the resulting 
equations in terms of expansion coefficients — {4ii\xj) 

i^j ~ £i) ^ij ^ y ' biaja + ^ ] {^ka^kiaj + £,kaCaikj} (18) 

a ka 

Here iijki is an anti-symmetrized two-particle matrix element iijki = tij^i ~ Ujik. Once the equations (^8|) are solved 
the "Brcit-dressed basis" can be determined as (jjj = (pj + '^[^ij't'i- The derived equations are essentially equivalent 
to the random-phase approximation or the self-consistent-field method, with an effective one-body Breit interaction 
serving as an external perturbation. The many-body diagrams for the amplitudes ^ij are shown in Fig. ^. By iterating 
these equations one sums a certain class of many-body diagrams to all orders in the Coulomb interaction. 

The resulting equations ( p^ are linear in the Breit interaction. It is worth noting that the Breit interaction, 
Eq. (^, is an approximation and terms nonlinear in the Breit interaction have no meaningful theoretical basis. 
Therefore linearized equations ( |l8| ) are conceptually more attractive than the self-consistent BCHF method based 
on an integration of Eq. (|p). However, the Breit contribution to atomic properties is relatively small and the much 
smaller terms nonlinear in the Breit interactions can be safely neglected at the present level of accuracy. 

An alternative approach to generating the BCHF basis set has proven to be more numerically robust and was 
employed in the present work. Two complete basis sets, CHF and BCHF sets, can be related by a unitary 
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FIG. 5: Separation of Breit correction into (a) two-body and (b) one-body (basis shift) parts. Here horizontal sohd and wavy 
hnes represent Coulomb and Breit interactions respectively; double horizontal line is a sum of Coulomb and Breit interactions. 
Vertical solid lines correspond to Coulomb-Hartree-Fock states and vertical double lines to Breit-CHF states. 



transformation 



Using Eq. (|T|) one determines expansion coefficients dkj and one-particle BCHF energies Sj from secular equations 

{sk - e,) dkj + Y,{^U)u d,j = , Vj . (19) 

i 

In this work the difference between the two Hartree-Fock potentials AJ7 was generated using finite-difference methods. 
The radial Coulomb-Hartree-Fock basis was approximated with B-splines pg|| and then transformed into the radial 



BCHF basis employing Eq. (19). Negative-energy states, £i < —rrieC^, were included in the diagonalization procedure. 

To summarize, third-order many-body calculations were performed in the Breit-Coulomb-Hartree-Fock basis with 
the two-body Breit interaction i?^^^ treated on equal footing with the residual Coulomb interaction. Sample many- 
body diagrams are presented in Fig. ^ and |l|. Contributions of negative-energy states, discussed for example in Ref. , 
were also included and found to be relatively small |]l0| . Two series of third-order calculations were performed, first 
with the Breit and Coulomb interactions fully included using the Breit-CHF basis set, and second in the CHF basis set 
without the Breit interaction and negative-energy states. The obtained differences are the Breit corrections analyzed 
in the following sections. 

Numerical calculations were performed using B-spline basis sets generated in a cavity of radius 75 a.u. This cavity 
size has been chosen for numerical consistency with the previous determination of parity-nonconserving amplitudes 
by Blundell et al. The numerical quasi-spectrum was represented by 100 negative- and 100 positive-energy states 
for each angular quantum number k. The intermediate-state summations were performed over 75 lowest-energy 
positive-energy states and 75 highest-energy negative-energy states for each partial wave Si/2 — ^ii/2- 



III. ATOMIC PROPERTIES 
A. Energies of valence states. 

At the Hartree-Fock level, the Breit interaction contributes less than 0.1% to all the energy levels considered in 
Fig. ||. Numerical values for Breit correction at the Hartree-Fock level are given in Table |; these were obtained as 
differences between one-particle energies in Breit-CHF and CHF approximations, i.e. — Sy. As in the traditional 
CHF calculations, the first-order many-body contributions to valence energies vanish identically in the BCHF basis. In 
the second order the corrections arise due to self-energy diagrams. For each valence state we perform two calculations 
with and without the Breit interaction and take a difference between the two values. Further, we distinguish between 
two classes of Breit modifications, one-body and two-body corrections, as illustrated for a diagram in Fig. |. The one- 
body contribution arises from a transformation of the one-particle basis from Coulomb-Hartree-Fock to Breit-CHF. 

The calculated values and breakdown on various contributions are presented in Table |. Apparently the Breit 
correction to the correlation part of the energy, 6Ebo, is equally important as the modification in the lowest order, 
(5i?HF- The interplay of Breit corrections to various many-body diagrams for the energy of the 75*1/2 state is remarkable. 
Here the two terms SEnp and SEbo,B'^^^ are almost equal and have opposite signs, the resulting modification being 
determined by relatively smaller two-body Breit correction SEbo, B'''^\ From Table | we see that generally the two- 
body Breit contributions are smaller than the one-body corrections; the -B*-^-* corrections become important when the 
cancellations are involved. 

The first study of correlated Breit corrections to the energies of Cs and other alkali-metal atoms has been performed 
in Ref. [^3l Based on formalism developed in Ref. the corrections have been determined as an expectation 
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TABLE I; Contributions of the Breit interaction to energies of valence electrons in cm^^. i?cHP are the energies in the Coulomb- 
Hartree-Fock approximation. ^-Ehf column lists corrections at the Hartree-Fock level defined as Ev — £v- Columns SEbo-, B'^^^ 
and SEbo, S'^' are the contributions in the second order due to one-body and two-body Breit interactions respectively. 



State Echf SE^f 5Ebo,B^^^ SEbo,B^^'> Total 

65^2 -27954 3^2 =^1:98 =1I83 ^2:6" 

751/2 -12112 1.1 -1.1 -0.28 -0.26 

6P1/2 -18791 7.5 -0.08 -0.28 7.1 

7Pi/2 -9223 2.7 -0.1 -0.10 2.5 

6P3/2 -18389 2.9 -1.8 -0.25 0.84 

7P3/2 -9079 1.0 -0.56 -0.09 0.38 

5PI3/2 -14138 -10.2 -12 -0.35 -22 

5PI5/2 -14163 -11.8 -14 -0.33 -26 



TABLE II: Comparison of Breit corrections to energies of valence electrons in cm ^ . 



State 



This work 



CCSD 1 




-1.1 




0.72 


6.9 




2.6 




0.29 


0.45 



Gaunt, no S'^' 



651/2 
75*1/2 
6P1/2 

7Pl/2 

6P3/2 



7P 



3/2 



-2.6 
-0.26 

7.1 

2.5 

0.84 

0.38 



-4 

9 
2 
2 




value of the Breit correction to the Coulomb-Hartree-Fock many-body Hamiltonian, Eq. (Jl 1|) , 

6E, ^ {^f,^\6BH'c\^^,^) . (20) 

Here I'l'f^) is the linearized coupled-cluster wavefunction limited to single and double excitations from the reference 
Slater determinant a|^|Ocoro) (CCSD method). The correlations are built into these wavefunctions. The CCSD formal- 
ism accounts for a complete third order of MBPT with certain classes of diagrams summed to all orders. However, 
the random-phase-approximation sequence of diagrams, important for the self-consistent treatment of the Breit inter- 
action, is missed starting from the fourth order. The approach employed by Kozlov et al. |^ is similar to the present 
method, but in Ref. |^ the full Breit interaction has been approximated by the Gaunt term and contribution due to 
effective two-body interaction B'^^^ has been neglected. 

Comparison of our results with the CCSD values and results of Ref. is presented in Table The best agreement is 
for states where there is no cancellation between corrections SE^p and 6Ebo, B'^^\ The values differ significantly 
for iS'1/2 states where strong cancellations, emphasizing higher orders of MBPT, are present. One should keep in mind 
that these discrepancies arise only when the Breit corrections, due to cancellation effects, are small and do not have an 



enhanced effect on evaluation of parity non-conserving amplitudes, discussed in Section IV. Calculations |^ are less 
consistent with the present and CCSD results, with discrepancies caused by approximation of full Breit interaction 
by Gaunt term and neglect of two-body part of the Breit interaction. 



B. Magnetic-dipole hyperfine structure constants 

To reiterate discussion in Section ||, we performed two series of computations: (i) traditional Coulomb and (ii) fully 
including Breit interaction. The difference between resulting values defines the Breit correction. From parametrical 
argument it is assumed that the leading Breit corrections arise from induced modifications of the dominant traditional 
Coulomb diagrams. Therefore the calculations were limited to the random-phase- approximation (RPA) and Brueckner 
diagrams (see Fig. |l|). Further the RPA sequence was truncated at the third order. 

In calculations of Breit corrections to ^■^'^Cs hyperfine-structure magnetic-dipole (HFS) constants A nucleus was 
modelled by a uniformly magnetized ball of radius Rm = 5.6748 fm. The gyromagnetic ratio for ^^'^Cs nucleus is 
gi = 0.73789 M. 



The breakdown of Breit corrections to various classes of many-body diagrams is given in Table III, Clearly the 
Breit correction to correlations ((5A^p^^\ <^^bo) equally important as the modifications in the lowest order (5^Hp. 
As an extreme case, almost entire Breit correction to HFS constants of 6S and 7S states comes from correlations. 
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TABLE III: Breit corrections to magnetic-dipole hyperfine structure constants A of ^ Cs in MHz. Column CHF lists Coulomb- 
Hartree-Fock values. Breit corrections to a class X of many-body diagrams of order A'' are designated as SA^- Total Breit 
correction ^Axotai is a sum of modifications 5A^. 



State CHF SA\^ SA™ SA^ ^ATotai 

6S1/2 1425.3 0.011 4.1 0.79 4.87 

75'i/2 391.6 -0.029 1.1 0.08 1.15 

6P1/2 160.9 -0.68 0.39 -0.24 -0.52 

7Pi/2 57.62 -0.23 0.14 -0.059 -0.15 

6P3/2 23.92 -0.06 0.10 -0.008 0.034 

7P3/2 8.642 -0.022 0.038 -0.0020 0.014 

5PI3/2 18.23 0.099 0.10 0.11 0.31 



TABLE IV: Comparison of contributions of Breit interaction to magnetic-dipole hyperfine structure constants A of ^'^ Cs in 
MHz. 





65'i/2 


7Si/2 


6P1/2 


7Pl/2 


This work" 


4.87 


1.15 


-0.52 


-0.15 


Kozlov et al. [|* 


5.0 


0.8 


-0.2 


0.0 


Sushkov Inl" 


4.6 


1.09 






Safronova et al. |33|]'^ 


-4.64 


-0.83 


-0.87 


-0.29 


Blundell et al. 'fjiq 


0.00 


-0.05 


-1.25 


-0.39 



" Third-order calculations in the BCHF basis, Ref. ||7|] . 
^ Pull Breit interaction is approximated by the Gaunt term. 
Analytical aZ expansion with Z = 55. 

Second-order calculations in the CHF basis. See Ref. for details. 
' RPA sequence of diagrams in one-body Breit interaction 



There is a cancellation of various contributions to the HFS constants of P1/2 and P3/2 states. For these states the 
contribution of higher-order diagrams not included in the present third-order analysis can become enhanced. At the 
same the total Breit corrections to Aqs and A75 are expected to be insensitive to higher-order contributions. 

A comparison of our results, partially published in Ref. 0], with other calculations is presented in Table IV. The 
correction to hyperfine constants is very sensitive to correlations: e.g., Ref. |37| found a numerically insignificant 
modification for Aqs, while Ref. |3^ determined the modification to be large (-4.64 MHz), and the approach 
reported here yields -1-4.87 MHz. In the calculation of Ref. |^ the correction was determined as a difference of 
the Breit-CHF and CHF values, however such approach misses two-body Breit corrections of comparable size. In 
Ref. jl^, Q a second order perturbation analysis was used for the Breit interaction, but the important relaxation 
effect discussed earlier was omitted. The present calculation incorporates all mentioned diagrams and is also extended 
to third order. Motivated by strong dependence of results 0, |3^ on many-body corrections, Sushkov derived 
an analytical expression for Breit correction to HFS constants of S states. His results for 65* and 75 states are in 
an excellent agreement with the present calculations. Kozlov et al. ^ used an approach similar to Ref. 0. There 
is a cancellation of various contributions for the P states; both higher-order diagrams and an approximation of the 
full Breit interaction by the Gaunt term in Ref. are the sources of discrepancies between our results for 6P1/2 and 
7Pi/2 states. 

Hyperfine structure constants sample atomic wavefunctions close to the nucleus and provide a unique way of testing 
atomic-structure calculations of parity-nonconserving amplitudes. In Table ^ we combine the Breit corrections with 
the results of ab initio all-order Coulomb-correlated calculations pTf and compare the results with experimental 
values. It is clear that the Breit corrections uniformly improve the agreement. In particular, the theoretical HFS 
constants are improved to 0.1% for 6S'i/2, 7SV/2j and lPi/2 states except for 6P1/2 where the discrepancy becomes 
0.5%. While the achieved agreement in Table [V| is encouraging, one should keep in mind omitted QED corrections and 
higher-order contributions in Coulomb interaction. For example, an estimate [|8[ |9j for hydrogen-like Cs ion results 
in a QED correction to HFS constant A of S'-states at a few 0.1%. Due to electron-electron interaction in atomic Cs 
QED corrections can be significantly modified. Correlated calculations of QED corrections would be beneficial for 
reaching 0.1% level of accuracy needed for interpretation of parity-nonconservation and also for understanding the 
role of high-order Coulomb diagrams at 0.1% precision. 
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TABLE V: Comparison of theoretical and experimental hyperfine constants A of Cs. All-order Coulomb-correlated values by 
Blundell et al. are supplemented with Breit corrections. Deviation from experimental values are placed in square brackets. 



Coulomb 


2291.00[- 


-0.3%) 


544.09[^ 


-0.3%] 


293.92[0.7%] 


94.60[0.3%] 


Breit 


4.87 




1.15 




-0.52 


-0.15 


Total 


2295.87[- 


-0.1%] 


545.24[- 


-0.1%] 


293.40[0.5%] 


94.45[0.1%] 


Experiment 


2298.16 




545.90 




291.93(2) 


94.35(4) 



C. Off-diagonal 6S — 7S hyperfine-structure matrix element 

Experiments 1^ on parity-nonconscrvation (PNC) in ^■^'^Cs determine the ratio of 65"— 75* PNC amplitude SpNC to 
vector transition polarizability /3. The value of /? is difficult to calculate reliably since it vanishes in the nonrelativistic 
limit. Following suggestion |Q, Bennett and Wieman |^ determined a supporting ratio of [3 to off-diagonal magnetic- 
dipole matrix element Mhf with a precision of 0.16%. Such an approach eliminates P from the analysis, but requires 
an accurate value for Mjjf. 

The quantity Mhf can be expressed in terms of the off-diagonal magnetic-dipole hyperfine structure constant Aqs-7S- 
This constant can be well approximated by a semiempirical geometric-mean formula ||4l| 



A'es-7S « VAesArs , (21) 

where Aqs and A75 are precise experimental hyperfine structure constants. The accuracy of this expression was 
investigated in Ref. Most recently, Dzuba and Flambaum |Q employed several many-body techniques 

of increasing accuracy in the Coulomb interaction between electrons and found that this geometric-mean formula is 
accurate to a fraction of 10~^. Here we extend their analysis and rigorously consider the additional effect of the 
correlated Breit interaction. The Breit correction to hyperfine-structure constants Aqs and A75 is in the order of 
0.2% and it can affect the sub-0.1% accuracy of the Coulomb analysis [Q. 

Breit corrections are relatively small. If Eq. (|2l|) holds, the following relation between the Breit corrections {6 A) 
has to be satisfied 



^A'6S-7s 1 f SAes , SAjs 



(22) 



^6S-75 2 t ^65 Ars 

As a result of the correlated Breit calculations we find SAqs-7S = 2.4 MHz. With Aqs-7S — 1120.1 MHz, the ratio 

5Aes-7S 



^"^^^ ''^ — 2.1 X 10 Using the Breit corrections to HFS constants Aqs and Ajs from Table [II, the semiempirical 



r.h.s. of the above equation is also 2.1 x 10"'^. Clearly the accuracy of the geometric-mean formula (|2l] ) is not affected 
by the Breit correction. Qualitatively this can be explained by a close proportionality of 6S and 7S wavefunctions in 
the vicinity of the nucleus, where the main contribution to the HFS constants of S states is accumulated. 



D. Electric-dipole transition amplitudes 



Calculated Breit corrections to reduced electric-dipole matrix elements of various transitions in Cs are presented 
in Table VI . We note that the Breit corrections to the random-phase approximation diagrams are small compared to 
the lowest-order and Brueckner-orbital corrections. Generally the total corrections are rather small (~ 0.1%), with an 
exception of 6Si 12 ~ 7Pi/2 electric-dipole matrix element. Using the ab initio all-order Coulomb-correlated value by 
Blundell et al. {'oSii2\\D\\lP^/2) = 0.279, and adding the Breit correction of 0.0019, one finds (6S'i/2| |7Pi/2) = 
0.281 in a better agreement with the 0.284(2) experimental value of Shabanova et al. The relatively large Breit 
correction is caused both by an accidentally small matrix element and by admixture into (6S'i/2|^|7/'i/2) from a 30 
times larger 7Si/2 — TPi/2 niatrix element. 
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TABLE VI: Breit corrections to reduced electric-dipole matrix elements for transitions between low-lying valence states of Cs 
atom. Column SD, % represents a ratio (in %) of the total Breit correction to a matrix element calculated in a third order of 
MBPT. See caption of Table |n| for description of other columns. 



Transition 


CHF 


<5-Dhp 


RPA 




(5 _D Total 


5D,% 


65*1/2 - 6P1/2 


5.278 


0.00035 


-0.00022 


-0.0011 


-0.00097 


-0.02 


65i/2 - 6P3/2 


7.426 


0.00078 


-0.00045 


-0.0014 


-0.0011 


-0.02 


65i/2 ~ 7Pi/2 


0.3717 


0.0018 


-0.00013 


0.00021 


0.0019 


0.5 


65i/2 — 7P3/2 


0.6947 


0.00059 


-0.00020 


0.00011 


0.00049 


0.07 


75i/2 - 6P1/2 


4.413 


0.0046 


-0.000067 


0.00038 


0.0049 


0.1 


75i/2 — 6P3/2 


6.671 


0.0019 


-0.000011 


-0.00030 


0.0016 


0.02 


75i/2 — 7P1/2 


11.01 


-0.0011 


-0.000050 


-0.0018 


-0.0029 


-0.03 


75i/2 — 7P3/2 


15.34 


0.0007 


-0.00013 


-0.0019 


-0.0013 


-0.009 


5D3/2 - 6P1/2 


8.978 


-0.0044 


-0.00035 


-0.0082 


-0.013 


-0.2 


5D3/2 - 6P3/2 


4.062 


-0.0028 


-0.00019 


-0.0035 


-0.0065 


-0.2 



IV. PARITY-NONCONSERVING AMPLITUDE 65i/2 75i/2 

The parity-nonconserving amplitude for the 6S1/2 transition in -'^'^^Cs can be represented as a sum over 

intermediate states mPi/2 

^ (75|i?|mPi/2)(mPi/2|i/w|65) 
^ y {7S\H^\mP,/2){mP,/2\D\6S) ^^^^ 

Here D and iJw are electric-dipole amplitudes and weak interaction matrix elements, and Ei are atomic energy 
levels. The PNC amplitude is expressed in units of 10~^^i|e|ao(— Qw/-^), where iV = 78 is the number of neutrons 
in the nucleus of ^^■^Cs and Qw is the weak charge. In these units the results of past calculations for ^^'^Cs are 
-E-PNC — —0.905, Ref. and E'pnc = —0.908, Ref. The former value includes a partial Breit contribution 
-1-0.002, and the latter includes none. The reference many-body Coulomb-correlated amplitude 

EpNC = -0.9075 (24) 

is determined as an average, with the partial Breit contribution removed from the value of Ref. The major 
difference between present and previous calculation of Breit correction to the PNC amplitude is an additional 
incorporation of effective two-body part of the Breit interaction and Breit corrections to the correlations. 

It is convenient to break the total Breit correction SEp^c into three distinct parts due to corrections in the weak 
interaction and dipole matrix elements, and energy denominators, respectively 

SE^^c = EpnciSH^) + EpMSD) + Ep^ciSE) . (25) 

For example, the modification of the PNC amplitude due to the Breit corrections to energies SEnS, SEmP^^^ '^^^ 
expressed as 

EpM^E) = '— [5Ees - 5E„,p,,,) 

+ c.c.(7S' ^ 65*) , (26) 

where the last term stands for the complex conjugate of the first term with 65 and 75* states interchanged. The 
Breit corrections to energies and dipole matrix elements were discussed in the preceding sections; here we focus on 
corrections to the weak matrix elements. 

The overwhelming contribution from parity-violating interactions arises from the Hamiltonian 

Gp 

Hw = —^QwPnuc{r)-f5 , (27) 
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TABLE VII: Breit corrections to matrix 
of Table for description of columns. 


elements of weak interaction for ^'^'^Cs 


in units of 10 ^^i|e ao( — ( 


Qw/N). See caption 


Transition 


CHF 




d(-f/w)R,pA 


(5(//w)bo 


5(-ffw) Total 


QSi/2 - 6P1/2 
6>S'i/2 — 7P1/2 
7Si/2 - 6P1/2 
TSii2 " 7P1/2 


0.03159 
0.01891 
0.01656 
0.00991 


-0.00010 
-0.000058 
-0.000053 
-0.000031 


-0.00015 
-0.000091 
-0.000081 
-0.000048 


-0.000029 
-0.000014 
-0.000013 
-0.0000061 


-0.00028 
-0.00016 
-0.00015 
-0.000085 



where Gf is the Fermi constant, 75 is the Dirac matrix, and pnuc{r) is the neutron density distribution. To be 
consistent with the previous calculations the Pnuc('') is taken to be a proton Fermi distribution employed in Ref. 
The slight difference between the neutron and proton distributions will be addressed in the conclusion of this section. 
The d omin ant contribution to the PNC-amplitude, Eq. (p3[), comes from intermediate states 6P1/2 and lPi/2- In 



Table VII we present calculated third-order Breit corrections to the relevant matrix elements of weak interaction. 
Apparently, the dominant part of the Breit correction arises from modifications at the Hartree-Fock level and in 
random-phase approximation (PRA). All the corrections add coherently, and we do not expect that omitted higher- 
order diagrams to be important. In fact, third-order RPA corrections are a few times smaller than those in the second 
order, hinting at a good convergence of the present technique. 

The Breit correction to the PNC amplitude is determined by combining in duced modifications in matrix elements 



and energy denominators. The required corrections are summarized in Table VIII. The tabulated dipole amplitudes 
are related to the reduced matrix elements in Table VI as 



red corrections are summarized in Table VIII 
ments in Table VI as 



Before proceeding to the correlated calculations, it is worth examining the Breit contribution to the PNC amplitude 
at the Hartree-Fock level. Most of the Breit contribution to the PNC amplitude can be determined by limiting the 
summation over intermediate states in Eq. ( |2^ ) to the two lowest valence P1/2 states: 6P1/2 and lPi/2- In the 
CHF approximation one then finds E'pnc — —0.6888 (90% of the total value). The lowest-order corrections to 
matrix elements and energy denominators calculated as differences between Brcit-CHF and CHF values are listed in 



Table VIII. The resultant BCHF corrections to i?pNC are: 



E^^ci^H^) = 0.0022 (0.32%), 
E%c{5D) = 0.0020 (0.29%), (28) 
E^l^iSE) = -0.0019 (-0.28%). 

Here the percentage values in parentheses are taken with respect to the DHF value of PNC amplitude. The sum 
of these three terms leads to — 0.0023. Inclusion of intermediate states beyond &P1/2 and 7Pi/2 leads to a 

small additional modification to SE-p^c of -0.00004. The obtained lowest-order result is in agreement with the 0.002 
correction found by BlundcU ct al. In addition to the lowest-order the Breit correction in Ref. |^ also contained 
small random-phase-approximation diagrams in the Coulomb interaction for matrix elements of the weak interaction. 
The two-body Breit interaction has been disregarded in Ref. In the following discussion we will include these 
omitted effects. Note that if experimental energies (which effectively include the Breit interaction) are used in the 
energy denominators of Eq. ([2^), then the i?pNc('5-E) term must be excluded and the total correction becomes twice 
as large: 5Ey>nc = 0.0042. 



With further examination of the modifications of individual uncorrelated matrix elements summarized in Table VIII , 
one notices the following. 

(i) Weak interaction matrix elements are each reduced in absolute value by 0.3%, which is directly reflected in a 0.3% 
correction to the PNC amplitude. 

(ii) Modification of dipole amplitudes is strongly nonuniform. There are substantial corrections only to the 
&S1/2 - lPi/2 (0.5%) and 7S'i/2 - QP1/2 (0.1%) matrix elements. The large 0.5% Breit correction to (65i/2|i:'|7Pi/2) 
provides partial resolution to a long-standing 1.5% discrepancy of spectroscopic experiment |45| and ah initio Coulomb- 
correlated calculations |7[ p6[ . 

(iii) The largest modification in the energy denominators is 0.1% for £"75 — Eqp] however, this leads to a 0.3% correc- 
tion -EpNC . As recently emphasized by Dzuba et al. , such large sensitivity of the resulting PNC amplitude to 
small variations in individual atomic properties entering Eq. (|2^) arises due to a cancellation of relatively large terms 
in the sum over states. 

It is well known that correlations caused by residual Coulomb interactions not included in the Hartree-Fock 
equations can lead to substantial modifications of the lowest-order values. For example, the weak matrix element 
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TABLE VIII: Breit corrections to electric-dipole amplitudes, weak interaction matrix elements, and energy intervals; SX, I = 
Xbchf — XcHF, and (5X, I + II + III are the differences in the third order of MBPT. Weak matrix elements are expressed in 
units of 10~^^i\e\ao{ — Qw/N) and energies and dipole amplitudes in atomic units. 





DOi/2 0-^1/2 


DOi/2 1-^1/2 


iJl/2 'J-ri/2 


71?, ,„ _ 7P, ,„ 
tJl/2 'J 1/2 


Hw, DHF 


0.03159 


0.01891 


0.01656 


0.009913 


5Hw, I 


-0.00010 


-0.00006 


-0.00005 


-0.000031 


SHw, I+II+III 


-0.00028 


-0.00016 


-0.00015 


-0.000085 


D, DHF 


2.1546 


0.15176 


1.8017 


4.4944 


5D, I 


0.0001 


0.00073 


0.0019 


-0.0004 


5D, I+II+III 


-0.0004 


0.00077 


0.0020 


-0.0012 


AE, DHF 


-0.041752 


-0.085347 


0.030429 


-0.013166 


5AE, I 


-0.000020 


0.000003 


-0.000030 


-0.000007 


5AE, I+II 


-0.000045 


-0.000023 


-0.000034 


-0.000012 



(65i/2|ffw|6Pi/2) is increased by a factor of 1.8 by correlations due to residual Coulomb interactions. As demonstrated 
in the previous sections, correlations are also i mpor tant for a proper description of the Breit corrections. Examination 



of the third-order corrections listed in Table VIIl reveals that the corrections to weak interaction matrix elements 
become three times larger than those in the lowest order. Using third-order matrix elements and second-order energies 
the following ab initio corrections are determined; 

EpNciSHw) = 0.0043, 
Ep^ciSD) = 0.0035, 
Ept^ciSE) = -0.0028. 

Thus the lowest-order corrections given in Eq. ( |28| ) are amplified in higher orders. The calculated Breit corrections 
to the PNC amplitude is expected to be insensitive to the omitted higher-order diagrams. Such uncertainty can 
arise if the calculated Breit corrections to leading classes of many-body diagrams cancel. Indeed, the calculated Breit 
corrections always add coherently for matrix elements of weak interaction. There are also no strong cancellations 
between various Breit corrections to the relevant dipole amplitudes {nSi/2\D\n' (see Table [V^ ). Corrections 
to energy denominators are also stable with respect to the omitted higher-order contributions. For example, in 



Section III A we found that the Breit correction to the energy of 7Si/2 valence state due to cancellations of calculated 
contributions is small (-0.26 cm^^) enhancing possible effect of smaller higher-order corrections. However, in the 
calculation of the term Ep]<;c{^E) this (unstable) correction substantially appears only in a combination with a 25 
times larger and stable Breit correction (7.1 cm^^) to the energy of 6P1/2 state. 

We further improve the accuracy of the calculation by combining all-order Coulomb-correlated matrix elements and 
experimental energy denominators tabulated in Ref. |^ with the present third-order Breit corrections. The results 
are: 

Epf^ciSHw) = 0.0047 (0.5%), 
Epf^ciSD) = 0.0037 (0.4%). (29) 
Ept^ciSE) = -0.0030 (-0.3%) . 

Here the values in parentheses are defined relative to the Coulomb-correlated PNC amplitude, Eq. (p^. The total 
Breit correction to the PNC amplitude is SEp^Q — 0.0054. This result was first reported in 0]. Similar correction of 
0.0053 was obtained by Dzuba et al. Kozlov et al. |^ found a 20% smaller correction of 0.004; this discrepancy is 
most likely due to the omission of the retarded part of the Breit interaction in the calculations ||9[ . 

If the experimental energies (which incorporate Breit corrections by definition) are used in the energy denominators 
the term i?pNc('5-E) should be excluded and the semiempirical Breit correction becomes SEp^'^' = 0.0084. A minor 
difference of our treatment of the Breit correction to -EpNC a-nd Ref. j^] is the interpretation of results of previous 
Coulomb-correlated calculations ^. Authors |j assert that the ab initio 0.6% correction ( 0.4% in Ref. j9j ) 
should be used to augment i^p^c, Eq. |2j. Our approach is to exclude term Ep-t^ci^E). The difference in the two 
interpretations arises due to the difficulty of accounting for higher order Coulomb diagrams. Some semiempirical 
"fitting" or "scaling" procedure is used in practice to mimic the effect of the omitted contributions. Since only 
energies are known with a very high precision from experiments, experimental energies play a central role in such 
analysis. For example, the experimental energies were employed in eight out of ten test cases in the scatter analysis 
in the Table IV of Ref. [| (Phys. Rev. D) based on Eq. M). 
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The above discussion demonstrates some arbitrariness encountered in the analysis of theoretical values and assign- 
ment of theoretical uncertainty when a semiempirical adjustment of ab initio values is attempted. We notice that 
authors of Ref. [|[ ^ argue that the theoretical uncertainty of the PNC amplitude is in the order of 1%, therefore a 
0.3% difference between the two different interpretations of the Breit correction is irrelevant at this level. We believe 
that the most convincing error estimate would be a repetition of the scatter analysis originally performed in Ref. 
but with a Breit correction to the involved quantities included. Such calculation is beyond the scope of the present 
work. However, we expect that the resulting uncertainty would be better than the original 1% assigned to the re- 
sults 1^, ^ because of the better theory-experiment agreement for dipole amplitudes (see Q ) and hyperfine-structure 
constants (see Table Additional QED and especially neutron skin effects can further modify the value of -EpNC 
in a way that cannot be mimicked by the suggested analysis. 

The parity violation in atoms is dominated by the Z-boson exchange between atomic electrons and neutrons. 
However the reference value Ep^Q, Eq. (p^), is based on the empirically deduced proton distribution. The difference 
between the proton and neutron distributions is visualized as the neutron "skin" or "halo" . Here we update our 
previous treatment of the neutron skin correction with the most recent data from the literature. This correction was 
estimated in Ref. j|] but was not included in the final value for the PNC amplitude. It can be shown that the neutron 
skin correction SEp^'^ does not depend on the electronic structure, therefore it can be parameterized as 

(30) 

Here Rp is the root-mean-square (rms) radius of proton distribution and A_R„p is the difference between rms radii 
of neutron and proton distributions. This expression can be easily derived from analysis by Fortson et al. [Q. 
From nonrelativistic nuclear-structure calculations Pollock and Welliver ji^ concluded ARnp/Rp = 0.016 or 0.022 
depending on the model of nuclear forces. The calculations ||5^, ^ of nuclear distributions were relativistic and 
the corrections as twice as large ARnp/Rp = 0.043 — 0.053 were found. Therefore nonrelativistic calculations led 
to SEp^,^/ Ep^c of -0.1% and relativistic determinations to -0.3 - -0.4%. The latter values are comparable to the 
experimental error bar of the PNC amplitude Q and unfortunately it is difficult to assess the accuracy of the 
nuclear-structure calculations. 

Here we propose an alternative analysis of the neutron skin correction allowing to estimate the error bar. Indeed 
Trzcinska et al. very recently deduced neutron density distributions from experiments with antiprotonic atoms 
and concluded that 

N- Z 
~A 

Here N and A are the neutron and the mass numbers. Although experimental data for ^'^^Cs do not enter the 
analysis ||l^ , a wide range of stable nuclei was investigated. Assuming that this relation holds for ^'^^Cs we find 
ARnp = 0.13(4) fm. For ^^^Cs Rp = 4.807 fm H leading to 

^-^^ = -0.0019(8) . (32) 

Therefore the neutron skin corrects the PNC amplitude by -0.2% with an error bar of 30%. This uncertainty con- 
tributes only 0.06% to an error budget of the observed weak charge, i.e. the proposed determination of the neutron 
skin correction will be adequate until the 0.1% level of overall accuracy is reached. 

Combining the calculated semiempirical 0.9% Breit correction with the reference Coulomb-correlated value, Eq.(^4|), 



ARnp ^ (-0.04 ± 0.03) -I- (1.01 ± 0.15) — fm . (31) 



and the neutron skin correction, Eq.(32), one obtains the parity- nonconserving amplitude 

£'pNc(^^^Cs) = -0.8974 X 10""i(-Qw/iV) ■ (33) 



In section pi C we concluded that the result for the off-diagonal hyperfine structure matrix element Mhf by Dzuba 



and Flambaum |44 is not affected by the Breit correction. Using their value of Mhf together with the experimental 
results |, m we arrive at 

^^PNc("^Cs) = -0.8354(33) x 10"" a.u. (34) 
From the values above, the observed weak charge is 

gw("^Cs) = -72.61(28)expt(73)theor. 

This value differs from the prediction [|| of the Standard Model Qff = -73.20(13) by 0.7cr, versus 2.5cr of Ref. [|, 
where a is calculated by taking experimental and theoretical uncertainties in quadrature. Here we assigned 1% 
uncertainty to the theoretical PNC amplitude, Eq. (^3). The deviation stands at 1.3cr if 0.4% theoretical uncertainty 
is assumed as discussed by Bennett and Wieman [Q^. Following Ref. 0|. similar conclusion has been reached in 
Ref. §1. 
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V. CONCLUSION 



In this work we presented a relativistic many-body formalism for treating correction from the Brcit interaction. 
Numerical evaluation of the Breit corrections to a number of properties of cesium atom were carried out. In particular 
we considered energies, hyperfine-structure constants, electric-dipole transition amplitudes and 6S—7S parity violating 
amplitude. We demonstrated that the Breit corrections to correlations are as important as the modifications at the 
lowest-order Dirac-Hartree-Fock level. This work supplements Ref. with additional numerical results. The present 
treatment has been based on third-order relativistic many-body perturbation theory. In a few cases we observed 
intricate cancellations between the lowest-order and higher-order corrections. These are the counterintuitive cases 
where the most advanced methods of many-body perturbation theory originally developed for residual Coulomb 
interaction will have to be employed to obtain an adequate description of a small Breit correction. 

Is it possible to test the accuracy of the theoretical treatment of the Breit contribution for alkali-metal structure? 
One could consider the nonrelativistically forbidden magnetic-dipole transitions nS 1/2 — n'Si/2- A second-order 
analysis demonstrated exceptionally large contributions from the Breit interaction and negative-energy states for 
such transitions; more accurate all-order calculations would be desirable. At the same time an accurate experimental 
value for the 75 — 65 transition in Cs is available Q . 

We determined Breit correction to parity non-conserving (PNC) amplitude of the 65 — 75 transition. The calculated 
correction resolves most of the discrepancy Q between the standard model prediction and atomic PNC determination 
of the ^^'^Cs weak charge. 

Breit correction to the PNC amplitude is one of the smaller contributions which needed to be addressed in order 
to reach the next level of accuracy in ab initio calculations required for refined interpretation of parity violation. 
In this work we also evaluated and constrained the neutron "skin" correction. As discussed in Ref. j8[ |9[ [ll), the 
remaining corrections which can contribute at a few 0.1% are due to higher-order many-body diagrams in the Coulomb 
interaction and QED corrections. 
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